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Read, April 17, 1820. 

WHATEVER difference of opinion there may be among astro- 
nomers, relative to the respective convenience of the different 
methods of obtaining approximate elements of the orbit of a Comet, 
there can be none, I think, relative to the superiority of M. Laplace's 
method of correcting the approximate elements when found. But 
this method may, as it appears to me, be considerably improved in the 
following manner. 
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The method as given in the " Mecanique Celeste/'* may some- 
times become impracticable, or at least, lead to tedious computations 
for want of a knowledge of the degree of exactness required in the 
approximate elements. M Delambre, after describing -f M. Laplace's 
method, remarks, " Cette methode est fort belle; si elle est peu longue, 
elle r6ussit presque toujours; elle n'est pourtant pas infaillible, et apres 
de tres longs calculs,on a ete force d'y renoncer pour la seconde comete 
de 1813/' 

Not being able to refer to the observations used for that Comet, it 
is not certain whence the difficulty arose, whether in obtaining the 
approximate elements or in the correction. It is most likely in the 
latter; and, if so, the method here proposed would have prevented the 
very long calculations, and probably would have removed the whole 
difficulty. 

M. Laplace's method of correction is as follows : 

He takes three observations of the Comet, separated from each 
other by considerable intervals ; and, with the approximate perihelion 
distance (p) and approximate time of perihelion, he computes three 
anomalies, (», »', »") and three radii vectores (r, r 7 , r") Let / — »=U, 
»" — v = U'. Let also /3, /3', /3" and v, *, v" be the respective heliocentric 
longitudes and latitudes. He computes these by means of the ob- 
served geocentric longitudes and latitudes, the Comet's distances from 
the sun, and the earth's distances from the sun. From the heliocentric 
longitudes and latitudes he computes the angles at the sun, subtended 
between the first and second places of the Comet, and between the 
first and third. These angles he calls V and V respectively. Then, 
if the approximate elements were correct, 

V=U and V'=U'. 

This not being so, he supposes m =■ U — V and m'=U' — V. 

* Tom. 1. p. 225 &c. + Delambre, Astr. Tom. 3. p. 386. 



191 



Then 1 . Changing only the approximate perihelion distance by a 
small quantity, he obtains new values of m and m\ 

2. Changing only the time of perihelion by a small quantity, he 
obtains other values of m and m'. 

From the relation of these values, he obtains two equations, by the 
solution of which the corrections of the perihelion distance, and of 
the time of perihelion are obtained. 

If these be not sufficiently correct, as will generally be the case, 
new corrections must be investigated. 

Every correction of the elements requires nine long operations, in 
which, to obtain the necessary exactness, it is requisite to use 
logarithms of seven places of decimals. 

By the alteration I am about to suggest, there will be required for 
the first corrections only three operations, in which seven places of 
figures will be necessary ,• in the other part of the process, five or even 
four places of logarithms will be sufficient ; and, in the subsequent cor- 
rections, only three operations to seven places of figures will be neces- 
sary. The repetition of the part of the process in which four or five 
places only were used, may in general be dispensed with. 

Instead of two hypothesis as in M . Laplace's method, by which two 
equations with two unknown quantities arc obtained, I obtain with- 
out any hypothesis, two equations in which the fluxion of the peri- 
helion distance and of the time from perihelion are the unknown 
quantities, in the followiug manner: — 

Let the true values of U and V be U + 3U and V + $V. 

Then U + 3U = V + W and U' +31/ = V + W. 

Now we can compute 8U, iV, &V, W so as to obtain two equa- 
tions of the form 

adp + bdt V— U 
a'dp + b'de=Y- U' 
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From which the. values of dp and dt maybe obtained, t being the 
time from perihelion. 

Let T — the time in the table of the Comet of 109 days when 
anomaly =v 

A = change of anomaly in that Table in one day at time T 

r 
* = T>f ... (2) 



Then 



d,,AdT~(bye q uat.2)g-^) 

by equal (1) log r+2 log cos \ » — log p 
hence d log r = sin 1" dv tan \ » 4- <JP_ 

P 
d logr , x __ r . dp 



or (V) 



d v tan ■§- v 



sin 1" a ? 

Let S, T, C, represent the Sun, Earth, and Comet 
respectively, and P the projection of the Comet on 
the plane of the earth's orbit. 

The processes to obtain the heliocentric longi- 
tudes and latitudes, are 

sin SCT = sinSTCxST , C&T = 180»— SCT-STC 
r 

sin CTP (geo. lat.) sin CST 
sin CSP (hel. lat.) = ^TSTC 

cos CST 

cos hel. lat. 
Comets hel. long. _ Earth's hel. long. + TSP 
Having thus obtained /S, /3' and ic, 1/ 



(«) 



(*) 
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According to Laplace's method of computation, using an auxiliary 
angle A, 

sin 2 A = cos* ^ (ft — /3) cos tcos «•' 
sin \ V = cos (i t + i t + A) cos (i «■ + i t— A) 
From these formulae, we may proceed to investigate dV 
log sin hel. lat. = log sin geo. lat. — log sin STC + log sin CST 
therefore, d log sin hel. lat. = d log sin CST - d CST cot CST 
But, log sin SCT = log sin STC + log ST— log r 
d log sin SCT = — d log r 

or d SCT = ' ~ tan SCT d log r = __ SCT 

sin 1" 
therefore d CST =. — d SCT = x tan SCT 

Hence d hel. lat. (*■) = x tan SCT cot CST tan * . . . • (c) 
d hel. long. (/3) = d TSP 
log cos TSP = log cos CST — log cos * 
hence d TSP = cot TSP tan CST d CST— cot TSP tan <rd*- 
or substituting from equation (c) and the equation preceding (c) 
d heLlong. (/S)=cot TSP tan SCT x (tan CST— cot CST tan 8 *)....(e) 
log sin A = log cos •§■ (fi' — (2) + \ log cos * + £ log cos i/ 
Hence 2dA=— tan A{d(3T— d/3)tan£(j3'--/3)+d*tan r + d* tan* (f) 
and dV=— \ (d«- + d*'+2dA) tan ± V tan (-L r + ±r' + A)) ( , 
— i (d* + d^— 2dA) tan £ V tan (£ t+ ^r' — A) J ^ 

Thus by equation (a) we obtain the value of d/ — d» = dU and by 
equations (a), (6) (c) (e), ^/J and (g), the value of dV and thence 
the equation U + dU=V+dV in terms containing only known t quan" 
tities and the two unknown quantities dp and dt. And in a similar 
manner the equation 

U' + dU' == V+ dV in terms containing only known quantities and 
the two unknown quantities dp and dt. 

The solution of these equations furnishes the values of dp and dt. 
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The coefficients of the unknown quantities in these equations would 
be correct, if we had the elements correct ; but, as the elements are 
only approximate, so these coefficients are only approximate, and their 
degrees of accuracy depend on the influence of the errors of the 
elements on these coefficients. 

Thus, if in computing dsr, SCT should be nearly a right angle, the 
corresponding error of da- would be great, compared to what it would 
be, if SCT were nearly 45° — In the former case, therefore, the 
elements would require to be much more correct than in the latter — 
If SCT should be so near a right angle, that it was uncertain whether 
it were actually greater or less than a right angle, that observation 
could not be used for the correction of the elements ; and so in similar 
cases. The magnitude of the angles SCT, CST, and TSP, will 
therefore point out t-he fitness of the observations used, and the degree 
of exactness required in the approximate elements. 

When Elliptic Elements are required, equations (a) and (b) are 
to be adapted to an elliptic orbit. Then another unknown quantity 
viz., dm, m being = axis major, 

will be introduced ; and so four observations at least must be used, 
as there will be three equations with three unknown quantities. 
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Application to the Comet of July 1819. 



The observations made on the north meridian, at the Observatory 
of Trinity College, Dublin, gave 



1819. 


Mean time of ob- 
servation at Obser- 
vatory T. CD. 


Longitude by 
observation. 


Latitude by 
observation. 

N. 


July 4 
5 
6 


12 h 7' 21", 8 
12 7 38, 4 
17 7 46, 4 


100 37 35, 7 

101 16 32, 3 
101 54 0, 


22« !4 53} 9 

23 33 1, 5 

24 39 30, 4 



The latitudes and longitudes are corrected* for aberration and 
parallax, and the longitudes are reckoned from the mean equinox. 

The method of Laplace-f- is very convenient for the first approxi- 
mation, when the observations are near together, and made with 
such exact meridian instruments as we possess at our Observatory. 

That method gave me, from the above longitudes and latitudes,^ 
perihelion dist. = ,3397, and time of perihelion June 27 d . 15\ 5'. 28". 



* A first approximation, which I computed merely to obtain the outlines of the orbit, enabled 
me to obtain the corrections for aberration and parallax. 

f Mec. Cel.tom. 1. p. 221, &c, 

£ The sun's places were computed from the Nautical Almanac, but the log. of the sun'* 
distance was taken from the Conn, des Temps, the log. of the distances in the Nautical Al- 
manac not being sufficiently exact to be used in the correction of the approximate elements. 
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By observations on the north meridian, July 4, 13, and 20, we 
have 



1819. 


Mean time of ob- 
servation at Obser- 
vatory T. C. D. 


Longitude by 
observation. 


Latitude by 
observation. 

N. 


July 4 
13 
20 


12 h 7' 21, 8 
12 4 29, 
11 55 13, 3 


100 37 35, 7 
105 41 6, 4 
108 45 50, 1 


22° 14 53, 9 
28 52 5, 3 
30 9 58, 5 



The observations of July 4 and 13, and the above approximate 
elements, give U— V=— 588",7. 

And those of July 4, and July 20, give U'— V'=— 655, 5. 

The processes for finding U, U' and V, V give the angles, by 
which the coeff. of dp and dt are determined. By help thereof we find 

i' » 

dv =18841 d«— 572072 dp 

d/ = 9380 At— 657490 dp 

d»"= 5859 dt— 595680 dp 



x = 



_ d log r __ 



sin 1 



7555 dH-377800 dp 



, d loff r' 
X = sin 1" = 7531 dt+79355 d P 



a!'= 



d log r" 



6065 dt— 9430 dp 

d/3=— 0,927* 
d/3'= + 11,886*' 



sin 1" 
dr = —4,682* 

dr = +0,429*': 

oV" = —0,641 2" : d/3" =+ 3,08 1 x" 
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By help of these quantities we find 

dV = 4,639 x+\,S6$x' 
dV'= 4,641^+1,036*" 
And finally we obtain 

58576 dt± 1986238 dp =—588,7 ") (p) 

54328 dt + 1 767 1 98 dp = —655,5 5 

And we deduce 

dp = 0,031463 
d/= —0, 05958 
Hence the first corrected perihelion distance = 0,341163 

And time of passage 

through perihelion, . June 27 d . 16 h 31' 16" 

With these quantities we again compute 

U— V =—39", 5 and U'—V = — 28",2 
Substituting these quantities on the right hand side of the equa- 
tions (p) ; we obtain again 

dp =—,0001 12 and dt = +, 003126. 
It is unnecessary to continue the approximations farther, and these 
give the perihelion distance and time of perihelion sufficiently exact. 

To find the place of the node, place of perihelion, &c. it is ne- 
cessary to know the hel. longitudes and latitudes according as 
they would result from the corrected element* now found. These 
are easily had by substituting in the above values of iW, d/3, &c. 
the second values of dt and dp, and applying the results to the re- 
lative longitudes and latitudes found in the second operation for find- 
ing V and V. The elements thus obtained from the observations of 
July 4, 13, and 20', are 
Passage of perihelion, 
Mean time at Observa- 
tory, Trinity College i 
Dublin » 
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June 27 d . 16\ 26'. 46" 
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Perihelion distance 0,341051 
Longitude of node 9'. 3". 43'. 44" 
Inclination 80. 45. 53 

Place of Perihelion 9.17. 5. 5 
Motion direct. 



Comparison of the observations with the Elements. 



1819. 


Mean time of ob- 
servation at Ob- 
servatory T. C D. 


Longitude by 
observation. 


Latitude by 
observation. 


Errors of elements 
in 


Long. 


Lat. 


July 4 
5 

6 


12 h 7' 21", 8 
12 7 38, 4 
12 7 46, 4 


100° 37' 35'', * 
101 16 32, 3 
101 54 0, 


22" 14' 53", 9 

23 33 1, 5 

24 39 30, 4 


-17", 2 
—15, 6 
—15, 3 


— 4", 
+ 4, 2 
—14, 8 


13 
14 
15 


12 4 29, 
12 3 28, 
12 2 21, 3 


105 41 6, 4 

106 9 27, 6 
106 37 10, 


28 52 5, 3 

29 9 37, 5 
29 24 26, 1 


+ 0, 6 
+ 5, 6 
— 1, 5 


— 0, 2 
-0, 7 

— 1, 1 


20 
21 

24 
28 


11 55 13, 3 
11 53 33, 2 
11 48 5, 4 
11 39 57, 9 


108 43 50, 1 

109 7 14, 4 

110 14 24, 1 

111 36 44, 9 


80 9 58, 5 
SO IS 7, 7 
30 25 52, 5 
30 33 3, 9 


+ 18, 9 
+ 20, 8 
+ 16, 9 
+20, 9 


+ 3, 3 

+ 8, 2 
+ 14, 9 
+13, 9 



